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Abstract 



Coulomb systems in which the particles interact through the <i-dimensional Coulomb 
potential but are confined in a flat manifold of dimension d— 1 are considered. The Coulomb 
potential is defined with some boundary condition involving a characteristic macroscopic 
distance W in the direction perpendicular to the manifold : either it is periodic of period 
W in that direction, or it vanishes on one ideal conductor wall parallel to the manifold at 
a distance W from it, or it vanishes on two parallel walls at a distance W from each other 
with the manifold equidistant from them. Under the assumptions that classical equilibrium 
statistical mechanics is applicable and that the system has the macroscopic properties of 
a conductor, it is shown that the suitably smoothed charge correlation function is univer- 
sal, and that the free energy and the grand potential have universal dependences on W 
(universal means independent of the microscopic detail). The cases d = 2 are discussed in 
detail, and the generic results are checked on an exactly solvable model. The case d = 3 
of a plane parallel to an ideal conductor is also explicitly worked out. 
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1. INTRODUCTION 

The present paper is about Coulomb systems, in which the particles interact through 
the (i-dimensional Coulomb potential (or some variant of it) but live in a space of di- 
mension d — 1 (this is what we call restricted dimension). Examples are particles in a 
plane interacting through the usual 1/r Coulomb potential, or particles on a line inte- 
racting through the two-dimensional — In r Coulomb potential. We are interested in the 
classical equilibrium statistical mechanics of some of these systems ( "classical" means that 
the quantum effects are disregarded), with the purpose of exhibiting universal properties 
("universal" means independent of the microscopic detail). This universality is closely 
related to the universality of the macroscopic electrostatics of conductors. Therefore, from 
the beginning, we assume that the systems under consideration have a macroscopic con- 
ducting behavior ; this excludes for instance a two-component one-dimensional log-gas at 
a too low temperature, since this gas is then insulating/ 1 ) 

The charge correlations have already been shown to be universal for several classes of 
conducting systems of restricted dimension. In the simple above-mentioned cases/ 2 ) the 
charge-charge correlation function S(r) is 



in a plane with 1/r interactions, and 



S(r) = S (1-2) 



7r 2 r 2 



in a log-gas on a line ; ks is Boltzmann's constant, T the temperature, and r the distance. 
(These expressions (1.1) and (1.2) are only macroscopically valid : the distance r must 
be large compared to the microscopic scale, and possible oscillations of the correlation 
function have to be smoothed away). Another occurrence of universality is for systems 
with a (i-dimensional Coulomb interaction confined in some appropriate (i-dimensional 
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domain ; then, not only are there universal surface charge correlations/ 2 - 1 but also the free 
energy and the grand potential have universal finite-size corrections. ( 3-5 ) Here we shall 
be concerned with infinite flat systems of restricted dimension d — 1, with some finite-size 
effect brought in by a boundary condition on the electric potential in the <ith dimension. 
Universal behavior will be found both for the charge correlations and the free energy or the 
grand potential. For d = 2, these generic properties will be checked on a solvable model. 

For instance, in Section 2, we consider a conducting log-gas on an infinite straight line 
(the x-axis) with the logarithmic interaction modified by the constraint that it be periodic 
of period W in the transverse direction y. Then, for macroscopic distances, the smoothed 
charge-charge correlation function is found to be 

= _ k B Tcosh(*x/W) 
W 2 smh 2 (7vx/W) 

while the free energy per unit legnth / and the grand potential per unit length u exhibit 
a finite-VF correction as W — > oo : 

f(W) = /(oo) + k B T^ + oiW- 1 ) (1.4a) 

oj(W) =u(oo) + k B T^- + o(W- 1 ) (1.46) 

o W 

These generic results (1.3) and (1.4) are checked in the special case of a one-component 
log-gas which is a (partially) exactly solvable model. 

Similar results are obtained for a conducting line at a distance W from an ideal 
conductor in Section 3, and for a conducting line with ideal conductors on each side of it 
in Section 4. An example of higher dimension is considered in Section 5 : a conducting 
plane at a distance W from an ideal conductor. 
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2. LOG-GAS ON A LINE WITH TRANSVERSE PERIODIC BOUNDARY 
CONDITION 



2.1. The system 

The two-dimensional Coulomb interaction between two point charges q and q' located 
in the xy plane at r = (x,y) and r' = (x',y') is qq'Go(r,r'), with Go(r, r') a solution of 
the Poisson equation 

AG (r,r') = -27r<y(r - r') (2.1) 



In the strip — W/2 <y,y'< W/2, a solution of (2.1) with periodic boundary conditions at 
y = -W/2 and y = W/2 is 



G (r,r') = -ln 



W . , tt(z-z') 
— sinh 

TV W 



(2.2) 



where z = x + iy and z' = x' + «?/ are the complex coordinates ; the factor VF/tt in 
(2.2) ensures that, in the limit W — > oo, one recovers the usual two-dimensional Coulomb 
interaction — ln|z — We consider some one-dimensional system of charges, on the x 
axis, with the interaction (2.2), where now z = x and z' = x' ; some short-range interaction 
might also be present. It may be noted that the interaction 



v(x — x') = — In 



— sinh — (x — x') 
ix W 



(2.3) 



on the x axis interpolates between a logarithmic interaction — In \x — x'\ at short distances 
and a linear interaction (one-dimensional Coulomb interaction) — (n/W)\x — x'\ at large 
distances. 

The system is assumed to have the properties which characterize a conductor, and to 
be globally neutral. The period W is macroscopic. 



2.2. Electric potential correlations 

For investigating the charge correlations and the thermodynamics, we first need infor- 
mation about the correlations of the electric potential. Let 0(r) be the microscopic electric 
potential created at r by the charges of the system and let us consider the correlation func- 
tion < (p(r)(j)(r') >. For r and r' at macroscopic distances from the x axis, this correlation 
function can be obtained by the method described in ref. 1, using linear response theory 
and the conducting behavior assumption, as follows. Let us put an infinitesimal test charge 
q at r'. Its interaction with the system is described by a Hamiltonian q(p(r'), and by linear 
response theory the average potential at r created by the charges of the system is changed 
by 

<ty(r) = -0q < 0(r)0(r') > T (2.4) 

where f3 = 1/ksT and < ••■ > T means a truncated statistical average (< AB > T 
=< AB > — < A >< B >) ; here however < 0(r) >= 0, and the mark "truncated" 
is superfluous. The total potential change at r is gG(r, r'), where G(r, r') is given by the 
macroscopic electrostatics of conductors, i.e. G is the solution of the Poisson equation for 
a point charge g at r', with the constraint that it is periodic of period W in y and that it 
vanishes on the x axis y = : 

sinh ^y(z-z') 



G(r,r') = - In 
G(r,r') = - In 



sinh 


2\y( Z 


-z>) 


cosh 


2\y( z 


-z>) 



cosh -z') 



if yy' > 
if yy' < 



(2.5a) 
(2.56) 



(G does not vanish for yy' < 0, i.e. the conducting line does not screen the regions 
y > and y < from each other ; this is an effect of the periodic boundary condition 
which connects these regions across the boundaries at y = ±W/2). That part of the total 
potential change which is created by the charges of the system is 



8cj ) (r) = q[G(r,r')-Go(r,r')} 



(2.6) 



From (2.4) and (2.6), one obtains for the correlation function 

(3 < 0(r)0(r') > T = q [G (r, r') - G(r, r')] 

in terms of (2.2) and (2.5). 



(2.7) 



2.3. Charge correlations 

Let a(x) be the charge per unit length on the x axis and let the macroscopically 
smoothed charge-charge correlation function be S(x-x') =< a(x) cr(x') > T . Since 2na(x) 
is equal to the discontinuity of the electric field component —d<J)/dy across the x axis, and 
since Go does not contribute to that discontinuity, one finds from (2.7) 



(3 S(x — x) = 2 lim 

y,y'^0 



<9 2 G(r, r') 



+ 



d 2 G(r, r') 



yy'>o 



dy dy' 



yy'<0 



(2.8) 



dy dy' 

Using (2.5) in (2.8) gives (1.3). 

An alternative derivation of (1.3) can be given by considering only the x-axis sys- 
tem, with the interaction (2.3) given, without any reference to the "outside world". One 
starts with the assumption that an external infinitesimal linear charge density qexp(ikx) 
is perfectly screened (for a wave-number k small enough, i.e. macroscopic). Therefore 
the system responds by creating a charge density — q exp(ikx). The interaction Hamil- 
tonian of the external charge with the system is qv(k) a(—k), where v(k) and &(—k) 
are the Fourier transforms of v(x) and a(x) respectively. By linear response theory, 
-q = -(3qv(k) < a(-k)a(k) > T , i.e. 



(3S(k) 



1 



v(k) 



(2.9) 



where S(k) =< a(—k)a(k) > T is the Fourier transform of S(x). Our definition of the 
Fourier transforms is, for instance, 

POO 

v(k) = dx e~ lkx v(x) (2.10) 
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Using (2.3) in (2.10) gives (in the sense of distributions) 



s 7T , Wk , j \ 

v(k) = -ctnh (2.11) 

k 2 



and, from (2.9), (1.3) follows. 

2.4. Free energy or grand potential 

The thermodynamic potential to be considered is the free energy per unit length / 
if the canonical ensemble is used, or the grand potential per unit length u> if the grand 
canonical ensemble is used. As a starting point we consider the derivative df/dW (at 
constant densities) or dui/dW (at constant fugacities). 

Let us draw some line parallel to the x axis ; this line divides the plane into two 
regions R on its right and L on its left. The derivative df/dW or du/dW is the force per 
unit length that region R exerts on region L, i.e. the T yy component of the Maxwell stress 
tensor : 



= T yy' = TZ< EyirY - E x (vY > (2.12) 



dW dW yy 4tt 

where E(r) = — V0(r) is the electric field at r (T yy should be independent of r). Since 
< E(r) >= 0, we can replace < • • • > by < • • • > T in (2.12), and using the derivative of 
(2.7) with respect to r and r' (with the limit r' = r taken at the end of the calculation) 
one finds 

e T »> = ^ < 2 - 13 > 

By integration of (2.12), one obtains (1.4). 

Alternatively, (1.4) can be derived by considering only the x-axis system with the 
interaction v(x). The partition or grand partition function depends on W through v(x), 
and deriving it with respect to W gives a statistical average of dv(x)/dW to be taken with 
the charge correlation function S(x) (since dv(x)/dW vanishes at x = 0, the self-part of 
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S(x) may or may not be kept) : 

Using (1.3) and (2.3), or (2.9) and (2.11), in (2.14), one recovers the result f3df/dW = 
l3du/dW = -n/8W 2 and (1.4). 

2.5. Solvable model 

The one-dimensional one-component plasma with the interaction (2.3) is a solvable 
model for special value(s) of the temperature. The model consists of particles of charge e, 
with a particle density rj, and a uniform background of charge density — erj which ensures 
overall neutrality. 

The correlations are known at the special temperature such that V := (3e 2 = 2 ; they 
are the same ones as for some quantum zero-temperature system. ( 6 ) From eqs. (3.13a) and 
(3.21a) 3 of ref. 6, 

2 



(3 S(x,x') = -2 



1 £(x) e(-x') -£{-x) £(x') 
W 6[ (0; e~ 2 ^ w ) 2 sinh(7r|x - x'\/W) 



+ 2r]5(x-x') (2.15) 



where 

£(x) = e-™ /2W 6>i Q + n V x; e~^ w ^ (2.16) 

(6*i is a Jacobian theta function). The charge-charge correlation function S(x, x') does not 
depend only on the distance \x — x'\ because the system has a crystalline structure. For q 
small, 

e x {u;q) ~ 2q 1 ^smu (2.17) 
For W large, (2.17) can be used in (2.15) and (2.16), giving 



3 A factor 2 is missing in the denominator of eq. (3.21a) in ref. 6. 
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P S(x,x') 



2 f n \ 2 ( n 

2 sin y— + Trrjxj sin y— — tvtjx 



'] e -ir(x-x')/W 



+2 sin 2 - nr]x) sin 2 Q + ttt/x') e *(*-*')/w - cos 2^ cos 2^'] ,x^x' (2.18) 

The macroscopically smoothed S(x — x') is obtained by averaging out the microscopic 
oscillations in (2.18), i.e. replacing the sin 2 terms by 1/2 and the cos terms by 0, which 
gives 



S(x - x') 



cosh[7r(x — x')/W] 
W 2 smh 2 [n(x-x')/W] 



(2.19) 



in agreement with (1.3). 



The free energy can be computed at the special temperatures such that T = 1, 2, 4, 
starting with the known partition function^ 7 ) for a finite system of N particles on a line 
of length L along the x axis, and an interaction which is periodic of period W in the y 
direction and also periodic of period L in the x direction ; on the line, this interaction is 



v (x — x') = — In 



0i 



7T(X - X') _ vW/L \ 

L ' 6 ) 



+ constant 



(2.20) 



(with a suitable choice of the constant, (2.20) goes to (2.3) as L — > oo). Including in the 
energy the particle-particle, particle-background and background-background interactions, 
one finds for the partition function 



Z N (T) = 



<(0;g) 



nt/2 



-N 2 r/8 



n a - « 2 ") 



,n=l 



-N 2 T/2 



m Cnt 



(2.21) 



where q := exp(— ttW/L) and CVr is the configuration integral 



N r L 



Cjvr : Yl / dxe I I 

e=i^° i<j<k<N 



(2.22) 
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Cat has been computed in ref. 7 for V = 1, 2, 4. In the thermodynamic limit iV — > 
oo, L — > oo, A^/L = 77 fixed, one finds for the free energy per unit length /(IT; W) = 
-limL-^nZjvfT) 



3f(l:W) = -r?ln^ + 



eft In 



22. e -27rW?7(n 2 +2nt) 



E 



(2.23a) 



0/(2; WO = -»7 ln(27r) + ^ - ry £ In (l - e " 4 ^) 



n=l 



(2.236) 



7T 



■jf 

^0 



j 9/(4;W) = -T / ln(87r^) + — -»7 / eft In 



In the large — W limit, in all cases, 



oo 

E 



(2n + t)e- 4lTW7l( - n2+nt) 



Pf(r;W)~/3f(r;oo) + 



TV 

8W 



pf(T;oc) = V 



1 - ^\ ln(2n V ) - ^ln^ + £ + In (I )! - ln(27r) - 1 



(2.23c) 



(2.24) 
(2.25) 



3. LOG-GAS ON A LINE PARALLEL TO AN IDEAL CONDUCTOR 



3.1. The system 



The region of interest is the half-plane y > 0. The electric potential is constrained to 
vanish on the X clXIS, 1. e. the x axis is an ideal conductor at zero potential. A solution of 
(2.1) with that boundary condition is 



G„(r,r') = -ln 



z — z 



z — Z' 



(3.1) 



We consider some one-dimensional system of charges, on the line y = W, with the corres- 
ponding interaction 

1 Cx-xM 2 

(3.2) 



v(x — x ) = — - In 



2 (i-x') 2 +4r 
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plus perhaps some short-range interaction. It may be noted that this interaction interpo- 
lates between — In \x — x'\ at short distances and 2W 2 / (x — x') 2 at large distances. W is 
a macroscopic distance, and the system is assumed to have the properties of a conductor. 

3.2. Correlations 

We follow the same steps as in Section 2. Let us first consider the case when the 
conducting system on the line y = W is kept at zero macroscopic potential. Now G is 
the solution of the Poisson equation for a point charge q at r' with the constraint that it 
vanishes on the lines y — and y = W : 



G(r,r') = - In 



sinh^(z -z') 
sinh M z -z') 



if < y, y' < W 



(3.3a) 



G(r, r') = - In 



z-z' - 2iW 



if y, y' > W 



(3.36) 



G(r, r) = if < y < W, y' > W ; or y > W, < y' < W 



(3.3c) 



From the analog of (2.8) one now obtains the universal correlation function 



p S(x) 



1 



1 



(3.4) 



2n 2 x 2 



8W 2 sinh 2 ^fL 



An alternative direct derivation of (3.4) uses (2.9) and the Fourier transform 




(3.5) 



3.3. Free energy or grand potential 



Still assuming that the line y = W is kept at zero macroscopic potential, we can use 
(2.7) (with < • • • > T =< • • • >), (3.1), and (3.3a) for computing the T yy component of the 
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stress tensor at some point r between the ideal conductor and the conducting line. One 
finds 

The VF-dependence of the free energy or grand potential is given by 

df du 



dW dW 



Tyy (3-7) 



provided / and uj are properly defined with a Hamiltonian which includes the self-energy 
interaction (l/2)g 2 ln(2W) of each particle of charge q with its image. It should be noted 
that, in the limit W — > oo, this self-energy and the two-body interaction (3.2) generate 
a well-defined total Hamiltonian containing only a two-body interaction — ln|x — x'\ ; 
therefore / and u are expected to have well-defined limits as W — > oo. From (3.6) and 
(3.7), one finds the large- W expansions 

Pf(W) = Pf(oo) ~^ + o {W- 1 ) (3.8a) 



with the universal finite- W correction — n/24W. 

Alternatively, (3.8) can be obtained from (2.14) by using either (3.2) and (3.4) (with 
the integral on x in (3.8) defined as its finite part), or (2.9) and (3.5). Since one must 
keep the self-contribution from the (l/2)ln[(a; — x') 2 + AW 2 ] part of (3.2), it is indeed 
appropriate to use in (2.14) the full correlation function S(x) or S(k) which includes the 
self part ; when understood as the Fourier transform of /3S(k) = l/v(k) in the sense of 
distributions, (3.4) does represent the full j3S{x). 
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3.4. Non-zero potential difference 

We now consider the more general case when the system on the line y = W is kept 
at a non-zero macroscopic potential $ ; correspondingly, there is on that line an average 
linear charge density 

$ 

a = —— (3.9) 

(3.9) is the equivalent of the familiar charge-potential relation in a plane condenser. 

Provided one defines the charge correlation function as Six) =< a(0)a(x) > T , where 
the T (truncated) sign is now relevant, the calculation of S(x) is unchanged and (3.4) is 
still valid. 

The free energy and grand potential however get additional terms. < E y (r) > is no 
longer zero, and (3.6) must be replaced by 

I3T„ := JL < Ey (vf - EJrf >= £ [< E„(r) 2 - E x (vf > T + < E„(r) >] 

From (df /dW) a = T yy one obtains 

l3 f(W) ~ «(») - ^ + l3n*HV = W (oo) - ^ + |?1 (3.11a) 

while from (9a;/9VF)<i) = T yy one obtains 

It is well known that the macroscopic electrostatic energy $ 2 /AuW = na 2 must come with 
different signs in / and in ui. 

3.5. Solvable model 



The one-dimensional one-component plasma on a line parallel to an ideal conductor 
is a solvable model/ 8 ' 9 ) at the special temperature V := f3e 2 = 2 ; it is a special case of 
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the more general models of two-dimensional plasmas with an ideal conductor wall^ 10 ) or 
between two ideal conductor walls. ( 5 ) One uses the grand canonical ensemble, with a fixed 
linear charge density — erj for the background and a fugacity £ which governs the particle 
density. The distance W between the system and the ideal conductor can have any value 
(not necessarily large) to start with. 

By a simple adaptation of the formalism in previous work/ 5 ' 8 ' 9 ) we obtain for the 
grand potential per unit length u> (including the background self-energy) 

(3u = -Tr ln(l + K) + 2-KrfW (3.12) 
where K is the continuous matrix 

4tvt]W 

K(x, x') = i( — — - (3.13) 

v ' s x-x' + 2iW v ' 

Using the Fourier transform which diagonalizes K 

r oo ( 2TvCe 4 ^ w - 2kW if k > 

K(k)= dx' e lk{x '- x) K(x,x') = I (3.14) 

J -°° 1 if k < 

and writing the trace (per unit length) as (2n)~ 1 J dk, one finds 

fa = - ^ l n [l + 2nQe A ^ w ~ 2kW ] + 2tv V 2 W (3.15) 

JO 2"7T 

For obtaining a large- W expansion (r/W ^> 1) of (3.15), we make the change of variable 
AurjW — 2kW = —u, split the integral on u into two integrals in the w-ranges (—AnrjW, 0) 
and (0, oo), and write 

1 r° 

2nr] 2 W = -—— dulne u (3.16) 

4nW J_ 4irvW 

This leads to the still exact expression 

1 



^=-7/ln(20- 



4nW 



[ du ln (l + -^—eA + [ du ln (l + 27rCe" u ) 

J-ATT-qW V 2yr C / JO 



(3.17) 
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Finally, when nW 3> 1, we can replace the lowest bound of the first integral in (3.17) by 
— oo and change u into —u, which gives, up to exponentially small terms, 

/3a; ~ -n ln(2<) - du In [ (l + JLe"") (l + 2^-") 

= -r/ ln(27rO — [ln(27rC)l 2 (3.18) 

The corresponding particle density is 

n = -C^(M~V+^M2<) (3-19) 

and the charge density of the system is e(n — n) = (e/AnW) ln(27r£). The system is neutral 
when ( = 1/2/T. Otherwise, its average electric potential is $ = e(n—i])W = (e/An) ln(27rC) 
and therefore, since f3e 2 = 2, 

£ = (3.20) 
2n 

($ contributes a term e$ to the chemical potential, as expected). Using (3.20) in (3.18), 
with /3e 2 = 2, shows that (3.11b) is verified in the present solvable model. (3.9) is also 
verified. Furthermore, for W — > oo, /3a; (oo) = — r/ln(27rC), n = n (independent of £), and 
the thermodynamic relation j3f = /3a; + (ha. On becomes /3/(oo) = — r/ln(27r), in agreement 
with (2.25). 

The formalism of previous work^ 5 ' 8 ' 9 ) also gives the correlation functions in terms of 
the continuous matrix g(x, x') defined in matrix notation as 

K 

9 = 3.21 

y 1 + K V ; 

The particle density is 

ra = 0(O,O) (3.22) 

and the charge correlation function is 

S(x) = -e 2 \g(x,0)\ 2 + e 2 n5(x) (3.23) 
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From (3.21) and (3.14), one finds 

, nN f°° dk lkx K(k) 



dk 



Akx 



In particular, 



o 2tv 1 + K(k) Jo 2tt 1 + e 2H/ ( fc - 2 ^) 



(3.24) 



(3.25) 



In the case of interest r]W ^> 1, (3.25) takes the form (3.19) and (3.24) can be rewritten 



as 



g(x,0) 



dk 



Akx 



27T 1 + e 2W(k-2^n) 



1 



1 



2nix 1 + e - 47TnW 



+ 



f 

Jo 



dk e 



ikx 



w 



2tt ix 2 cosh 2 W(k-2>irn) 



(3.26) 



where an integration by parts has been performed. For nW ^> 1, up to exponentially small 
terms, by extending the last integral in (3.26) to the range (— oo, oo), one finds 



1 

■- + 



^2'Kinx 



x ■ ^L s inh^ 

7T 2W 



(3.27) 



The microscopic correlation function S{x) is obtained by using (3.19) and (3.27) in (3.23). 
It has oscillations of period n _1 , the average interparticle distance. When these oscillations 
are averaged out, (3.23) with fie 2 — 2 agrees with the universal form (3.4). 

4. LOG-GAS ON A LINE BETWEEN TWO IDEAL CONDUCTORS 

4.1. The system 

The region of interest is the plane strip < y < W. The electric potential is cons- 
trained to vanish on the lines y = and y = W, i.e. these lines are ideal conductors at 
zero potential. A solution of (2.1) with that boundary condition is 



G (r,r') = -ln 



smh^(z-z') 



(4.1) 
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We consider some one-dimensional system of charges, on the line y = W/2, with the 
corresponding interaction 



v(x — x') = — In 



TV 



taiah — (x-x') 



(4.2) 



plus perhaps some short-range interaction. Now v interpolates between — \n\x — x'\ and 
2exp[— (tt/W)\x — x'\]. Again W is macroscopic and the system is assumed to be a con- 
ductor ; it is kept at some potential <E>. 

4.2. Correlations 



Now, using 



G(r,r') = - In 



sinh^z -z>) 



W 

if < y, y' < — 



(4.3a) 



WW w w 

G(r,r') = if <y< —, — <y > <W ; or — <y<W, 0<y' < — (4.36) 

and using the analog of (2.8), one finds the universal correlation function 

pS{x) = -^-5 (4.4) 

W 2 sinh 2 ^ 

Alternatively, one can use (2.9) and the Fourier transform 

/ 7 \ 7T , Wk , A „. 

v(k) = T tanh — - (4.5) 
k 2 



4.3. Free energy or grand potential 



Now the stress tensor component T yy is found to be such that 



^ =8^ + l£ <4 ' 6) 
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The charge density is 



2$ 



Thus, 



and 



N ~ WW - — + ^j- = 0/(oo) - — + ^ (4.8a) 



4.4. Solvable model 

The one-component plasma is again a solvable model, with the present boundary 
conditions, when V := (3e 2 = 2. The formalism of ref. 5 is applicable, with now 

A 1T V W 

K(x,x') = ^- ^— -y (4.9) 

V ; 2W cosh[^(x-x')] K 1 



and 



K{k) = tt( - — (4.10) 
coshWk 



Thus 



r r dk i 



pTTTjW 

1 + nC-' 



+ \nrfW (4.11) 



cosh VF/c 

For obtaining a large-VF expansion of (4.11), we rewrite it as 

(3iv = - —In _ ^ +/ — lne^^ 4.12 

Jo 7T 1 + e IWk Jo n 

extract from the first integral in (4.12) the factor 

io ^ 1 ^ 24W 1 ' 

and regroup the other terms into two integrals on the /c-ranges (0, nrf) and (ttt}, oo). After 
simple changes of variable, neglecting exponentially small terms (which allows to extend 
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one of the integration ranges) , one obtains 

i f°° r / i 

~ ln(2. C ) + — - — jf A, In [(l + 2* e -) (l + — e- 
= -"'"(2-0 " ^7 " H2<)f (4.14) 

The corresponding particle density is 

n = -C-^(Poj) ~ 47 + -i- ln(27rC) (4.15) 

Neutrality still occurs when 2tc( = 1, (3.20) is still valid, and (4.8b) is verified by (4.14). 

The correlation function is still given by (3.21) and (3.23) now with (4.10). Thus 

r°° dk f°° dk cos kx 

9{x, ) ~ ^ l + J~ g 7rr]W cosh Wjt ~ y o Vl + 2e- 7r " w/ coshV17/ ; ; 

sin(7mx) 



Wsinh^ 



(4.16) 



and 



_l-cos(27m^) 
V ' W 2 sinh 2 ^ V ' ' 

When the oscillations of (4.17) are averaged out, the universal expression (4.4) is verified. 

5. COULOMB SYSTEM IN A PLANE PARALLEL TO AN IDEAL 
CONDUCTOR 

5.1. The system 

The universal properties discussed in Section 2 to 4 can be generalized to systems of 
higher dimension. As an example, in three-dimensional space xyz, we consider a conducting 
classical system of charges confined in the plane z = W and kept at an average potential 
$, while the plane z = is an ideal conductor at potential zero. Such a model might 
be relevant for describing electrons trapped at the surface of liquid helium in front of 
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an electrode located under the surface. ( n ) We transpose the derivations and results of 
Section 3 to the present system. 

We note a position as r = (x, y, z) = (p, z) where p = (x, y). In the half-space z > 0, 
the Green function G 0: solution of 

AG (r,r') = -47r<5(r - r') (5.1) 
constrained to vanish on the plane z = and at infinity is 

Go(r ' r ' ) = ^-^1 (5 - 2) 
where r'* = (x ; , y' — z') is the image of r'. In the plane z = W, the interaction is 

v(p > p>) = \p~^7\ ~ W^WTWW (5 ' 3) 

5.2. Correlations 

The Green function which vanishes on both the z = and z = W planes is 

oo 

G(r,r')= G (r,r' + 2nWu) i£0<z,z'<W (5.4a) 

n= — oo 

where u is the unit vector along the z axis, 

G(r,r') = — ?— — if z,z'>W (5.46) 

|r-r'| [(p- (9 ')2 + ( z + z '-2W) 2 ] 1/2 

G(r, r') = if < z < W, z' > W ; or z > W, O < z' < W (5.4c) 

The potential correlation function is given by (2.7), and from the analog of (2.8) one finds 
for the surface charge cr(r) correlation function in the plane z = W 

47r n =o KP - P ) + (2nW) 2 ] x 
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At short distances, (3S ~ — l/47r 2 (p — p') 3 , in agreement with the formula (1.1) for a 
plane alone ; at long distances, the Euler-MacLaurin summation formula gives j3S ~ 

An alternative direct derivation of (5.5) uses (2.9) and the two-dimensional Fourier 
transform 

= jjVl - e -^M) (5.6) 



5.3. Free energy or grand potential 

In the region < z < W, using (2.7) gives for the T zz component of the Maxwell 
stress tensor, the force per unit area, 

where ((3) = 1.202 ... is a value of the Rieman zeta function. The last term of (5.7) 
is the standard attractive force between the plates of a plane condenser. The (usually 
much smaller) thermal term kBT((3)/8TrW 3 can be obtained as the classical limit of the 
celebrated more general theory^ 12 ) of Van der Waals type forces between macroscopic 
bodies. 4 

By integrating (df/dW) a = (duj/dW)$> = T yy , one obtains 

MW) = moo) - ^ - ^ (5.86) 

Alternatively the thermal term of (5.8) can be obtained by using (2.9), the analog of 
(2.14), and (5.6). 



4 Our result agrees with eq. (5.5) of ref. 12, with eq = oo and the integral evaluated 
exactly. 
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For the <i-dimensional analog of the present system, the thermal part of the force per 
unit area is found to be 



The occurrence of universal properties in conducting classical Coulomb systems is 
especially visible in systems of restricted dimensionality. The simplest, already known, 
example is the smoothed charge-charge correlation function S(r) for particles in a plane 
with 1/r interactions, as given by eq. (1.1). If (1.1) is understood in the sense of distribu- 
tions, the prescription for regularizing its integral is 



and therefore (1.1) even grossly represents the full S(r), including its self-part term, since 
(6.1) correctly expresses the screening rule. 

In the present paper, we have derived universal smoothed charge correlation functions 
for other geometries involving some boundary conditions for the electric potential outside 
the system of charges. These boundary conditions involve some macroscopic length scale 
W, and we have exhibited universal W- dependences of the free energy and the grand 
potential. 

In Sections 2 to 4, we have considered one-dimensional systems with a two-dimensional 
Coulomb interaction, because, for these systems, we had at hand exactly solvable models 
on which we could test our generic results. However, similar generic results can be easily 
derived for any dimension, and in Section 5 we have given an example involving the usual 
three-dimensional Coulomb law. 

A by-product of Section 4 is a new one-dimensional solvable model : particles on a line, 
interacting through the potential — e 2 In | tanh[(7r/2W)(x — x')]\, without any background 



k B T(d-l)T(l) ((d) 
ix d / 2 (2W) d 



(5.9) 



6. CONCLUSION 




(6.1) 
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nor self-energy. When T := f3e 2 = 2, the correlations and the thermodynamics can be 
obtained exactly, by a minor adaptation of Section 4. 
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